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Abstract 



Several nonlinear stochastic differential equations have been proposed in 
connection with self-organized critical phenomena. Due to the threshold con- 
dition involved in its dynamic evolution an infinite number of nonlinearities 
arises in a hydrodynamic description. We study two models with different 
noise correlations which make all the nonlinear contribution to be equally 
relevant below the upper critical dimension. The asymptotic values of the 
critical exponents are estimated from a systematic expansion in the number 
of coupling constants by means of the dynamic renormalization group. 
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I. INTRODUCTION 



Recently a lot of attention has been paid to the phenomenon known as self-organized 
criticality (SOC). Bak, Tang, and Wiesenfeld (BTW) Jl[ studied a cellular automaton model 
as a paradigm for the explanation of two widely occurring phenomena in nature: 1/f noise 
and fractal structures. Both have in common a lack of characteristic scales. This scale 
invariance suggests that the system is critical in analogy with classical equilibrium critical 
phenomena, but in self-organized criticality one deals with dynamical nonequilibrium statis- 
tical properties. On the other hand the system evolves naturally to the critical state without 
any tuning of external parameters. 

Several cellular automata models exhibiting self-organized criticality have been reported 
in the literature. In the original sandpile model of Bak. et al. [|l| the system is perturbed 
externally by a random addition of sand grains, once the slope between two contiguous cells 
have reached a threshold condition sand is transferred to its neighbors in a fixed amount. 
Taking this model as a reference, different dynamical rules have been investigated leading 
to a wide variety of universality classes: directed flows @, threshold condition imposed on 
the height, on the gradient or even on the laplacian ||, continuous variables with a full 
transfer from a cell instead of a fixed discrete amount deterministic perturbations in 

a nonconservative system anisotropy as some examples. 

On the other hand some authors have attempted to connect these cellular automata 
models showing SOC to nonlinear stochastic differential equations |T2]-[l^|. These continu- 
ous descriptions are built according to the symmetry rules obeyed by the discrete models. 
Anomalous diffusion equations with singularities in the diffusion coefficient have been con- 
sidered in order to study the dynamics of the avalanches generated in the critical state 
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In this paper we study two nonlinear stochastic differential equations whose deterministic 
parts keep the symmetries of the usual sandpile BTW models |7|JT7||. However, they have 
different noise correlations from which one gets different critical behavior. Both models 
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are treated independently and the noise correlations are obtained by different limits of a 
general Ornstein-Uhlenbeck process and in some way they represent different aspects of 
SOC. For both models an infinite number of nonlinearities is relevant below their upper 
critical dimension d c = 4, 2; by taking a finite number of nonlinearities we estimate the 
values of the critical exponents from an extrapolation. For one of the models we get the 
dynamical exponent obtained in the numerical simulations of systems exhibiting SOC. 



II. STOCHASTIC DIFFERENTIAL EQUATIONS 

The nonlinear stochastic differential equation we are going to study is 

= D ^E(r,t) + £ \ 2n+1 V 2 E^m + r](r,t) (1) 

01 n=l 

where D is the bare diffusion coefficient and \ are the coupling constants accounting for the 
nonlinearities of the model. The variable E which we call energy can have different physical 
interpretations || . The difference among both models lies in the noise correlation 

r 

< r}(r,t)r)(r',t') >= 2-5 d (r-r'). (2a) 

< r){r,t)r)(P,t') >= 2FS(t - t')5 d {f-r'). (2b) 

T being the intensity of the noise and T is a macroscopic time to be discussed in the last 
section. 

From Eq. ([!]) we can see that both models are translational and rotational invariant 
and that the energy is conserved in its deterministic part. An additional symmetry (E — > 



— E invariance) has been introduced since in [jrj] it is showed that it is generated by a 
renormalization group procedure when Eq. ([]]) is derived from an equation involving a 
threshold condition. 
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III. DYNAMIC RENORMALIZATION GROUP 



The scaling behavior of the system described by the nonlinear stochastic differential 
equation is obtained from the correlation function 



< (E(f , t ) - E(f + r, t + t)f r*f{t/r z ) 



(3) 



where x an d z are the 'roughening exponent' and the dynamical exponent respectively. The 
relevance of the different coupling constants can be checked by naive dimensional analysis: 
a change of scale r — > e l r is followed by t — > e zl t and E — > e xl E in order to (^) be satisfied. 
Under this scaling transformation z and x are chosen to keep the linear equation invariant. 
With these values one can see that all coupling constants are of the same order in e = d c — d 
and hence all nonlinear terms are equally relevant for d < d c . The upper critical dimension 
d c are 4 and 2 for models (pa]) and (|2b"|) , respectively. 

The dynamic-renormalization-group procedure consists of two steps: an elimination of 
the faster short wave-length modes followed by a rescaling of the remaining modes flTB,T§ 



The recursion relations after an infinitesimal transformation for the coefficients are [|17 
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(5b) 
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for models (pa]) and (pbj), respectively. In the previous equations 7j are numerical factors 
and bij are natural numbers accounting for the vertices Aj entering the renormalization of 
the vertex A n in the j-th term of the sum. 

Since all terms in this expansion are equally relevant below the respective upper critical 
dimension a calculation of the critical exponents would involve all the contributions. Never- 
theless the exponents can be obtained to first order in e = d c — d as a series in the number 
of coupling constants which are taken into account. For details of the calculations for model 



see 



|l7f where one obtains 

z = 2 - 0.32(4 - d); X = 0-18(4 - d) for d < d c = 4 
z = 2; X = ¥ for d > rf c = 4 



(6a) 



On the other hand the results for model fl2"15|) can be obtained in an equivalent way to be 
z = 2- 0.33(2 - d); x = 0.34(2 - d) for d < d c = 2 



z = 2; X = ^ for d > d r = 2 



(6b) 
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IV. RELATION WITH SELF-ORGANIZED CRITIC ALITY 

Up to now it is not clear yet the relation between the cellular automata models exhibiting 
self-organized criticality and the nonlinear stochastic differential equations analyzed in the 
present work. In the original sand-pile model, and in others derived from it, the system 
is perturbed externally until a site reaches a threshold condition and then an avalanche is 
triggered. During the evolution of the avalanche there is no external perturbation onto the 
system. This introduces two time scales in the model: a fast time scale (the diffusion scale) 
and a slow one (the noise scale). 

In order to be suitable for an analytical treatment one can assume the external noise to 
be described by a Gaussian process with zero mean and correlation function given by 

2F 

< r)(f,t)r}(r',t') >= — e'^'^S^f - r>). (7) 
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This is an Ornstein-Uhlenbeck process, with r the correlation time. This general case can 
describe two different limits depending on whether r is a microscopic or a macroscopic time. 
In the former case we get ( |2b"|) whereas in the later one obtains fl2a|) with a intensity of the 
noise that scales with 1/r. 

When analyzing the SOC models within a mesoscopic time scale we deal with the dy- 
namics of single avalanches. In this case model (|2a| ) is more appropiate since the noise scale 
is larger than the time scale we are interested in. A quenched noise of low intensity does 
not affect the evolution of the avalanche since it only depends on whether a site is above 



its threshold. This is the reason for which the estimation of the dynamical exponent (pa]) 
agrees with the values obtained in the numerical simulations and from scaling arguments 0] 
concerning the evolution of single avalanches. Nevertheless it can modify the roughness of 
the interface since it can make two distant sites to be correlated in an artificial way; however, 
the effect of the nonlinearities is to lower the value of the roughening exponent with respect 
to the linear model. 

SOC models can be analyzed from a different point of view by considering the avalanches 
to be instantaneous. In this case one is mainly interested in the form of the interface after 
the relaxation events (avalanches) have taken place and not in its dynamical evolution. 
This situation is properly described by Eqs. (|l|) and ( 2bD but the parameters entering the 



deterministic evolution diverge since they are related to the diffusion coefficient. Only finite 
values of the parameters make the problem suitable of an analytical treatment; on doing so 
we make the two time scales discussed above to be of the same order |fL4}| . Although it is not 



the situation one deals with in the simulations it is interesting to notice that the roughening 
exponent x is lowered with respect to the previously discussed model. Actually, one should 
expect a roughening exponent x < due to the fact that the energy at each site is bounded 
after the avalanches are finished. But as the two dynamics (diffusion and noise) overlap we 
can still get a rough interface, at least for low dimensionality. 

To summarize, we have analyzed two nonlinear stochastic differential equations by means 
of the dynamic renormalizaton group. Both equations share the deterministic evolution but 
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differ in the noise correlation. Due to the symmetries of the cellular automata models 
they are derived from, an infinite number of relevant nonlinearities must be taken into 
account. By taking a finite number of these nonlinear couplings we extrapolate to infinite 
and get an estimation of the critical exponents. These two models permit to describe 
different, and complementary, characteristics of cellular automata models exhibiting SOC. 
The quenched noise model reproduces the dynamical exponent obtained in the evolution of 
single avalanches but introduces additional correlations in the interface profile. On the other 
hand, the delta-correlated noise model gives a roughening exponent which is closer to what 
one would expect in the profile of the interface after the avalanches have taken place. 
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